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The main goal of this investigation was to apply a model of the 
reversing geodynamo, near the Cretaceous-Tertiary boundary. 

The model is based on the assumptions of non axialsymmetric dipole 
sources within the core. The Cretaceous-Tertiarv boundary (KTB) was 
sampled and magnetostratigraphically analysed, by R. Rocchia et al. 
(1990), Lowrie and Alvarez 1977 (Gubbio-section, Italy) as well as by 
Zeissl and Mauritsch 1989 in Gams (Austria). The available VGP-paths 
show only little differences in longitude when passing the equator and the 
transition can be called nearsided. Assuming that the behaviour of the 
geomagnetic field during a polarity transition is dominated by non-dipol 
terms (axialsymmetric and not axialsymmetric ones), we applied a 
flooding model to the observed data (similar to that of K. A. Hoffman). 
For octupole dominated transition-fields we start the onset of the reversal 
at high latitudes in the core, for quadrupole transition-fields in the 
northern hemisphere and the Gauss coefficients will be calculated for each 
field configuration. 


The Geological frame 

Since the magnetostratigraphic results of Gubbio (Italy) are already 
published (Lowrie and Alvarez 1977, Rocchia et al. 1990) in the present 
paper just information from the locality Gams are given. 

The Gosauian basin of Gams belongs to the eastern part of the 
Northern Calcareous Alps in Austria. From a local tectonic point of view 
it is part of the upper austroalpine Otscher-Decke and can be divided in 
a western and an eastern part. The eastern part, were the magneto¬ 
stratigraphic profil in discussion belongs to, is surrounded by an early 
Triassic “high” in the west, jurassic Plassenkalk in the north. In the South 
and East the boundary of the basin is formed by an overthrusting zone 
(H. Kollmann 1964 a, b) (Fig. 1). 

The gosauian sediments start in the eastern part of the basin of Gams 
with upper Santonian overlayn by the so called "deeper marl complex” of 
lower to early Upper Campanian age. This complex is followed by the 
Nierntaler beds of late Upper Campanian to Lower Paleocene age, the 
breccia-sandstone complex of Middle Paleocene and the clay-marl 
complex of Upper Paleocene. 
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Caused by the Lower Campanian overstrusting of the southern rim, 
Gosauian sediments were cut of in the Upper Campanian conglomerate 
and marl horizon and slightly transported to the north. 

The studied profile belongs to the Nierntaler beds, which consist of 
limestones and marls particularely in the Maastrichtian, marls, sand- and 
siltstones in the Danian. The whole flyschoid sequence is cut by several 
minor faults with displacements of 5-50 cm but always absolutely 
reconstructable. Sedimentological investigations (Preisinger 1990) 



Fig. 2: Scatch of the transition sequence and sample position. 
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described the sediment as being partly of turbiditic and partly hemipelagic 
origin. 


Sampling and rockmagnetic investigations 

The sampling was carried out by a portable drill in the Maastrichtian 
and in the sand/siltstone layers of the Danian. The weak to uncon¬ 
solidated parts of the Danian were sampled by a special tool, which allows 
to push a oriented plastic tube into the soil. The tube is marked in the usual 
paleomagnetic manner with an arrow and a fiducial line and is closed by 
a lid to prevent drying out of the sample. The extent of the disturbing 
effects of this sampling technique depends on the grain size of the sample 
material. The sampling interval varied between 5 and 15 cm for the whole 
K/T-boundary section, whereas for the detailed polarity section 1.0- 
1.5 cm (Fig. 2) were chosen. The polarity boundary from G“ to was 
chosen since the facies, the possibility of drilling and the lack of faults 
were good pre-conditions. 

On a large number of pilots, the rockmagnetic behaviour was 
studied. As it can be seen in Fig. 3 magnetite as well as haematite and 
maghaemite can be identified. Up to 150 °C to 200 °C there is also some 
influence be ironhydroxide which particularly influences the IRM- 
acquisition in the natural stage. After heating the pilots to 200 °C and 
repeating the IRM-acquisition there is always a remarkable shift towards 
soft, magnetite like behaviour (Fig. 3 e). 

Anblocking temperatures at about 650 °C and saturation at about 
0.5 T is a strong argument for maghaemite whereas hematite is 
characterized by unblocking temperatures of about 650 °C and saturation 
above 1.5 T. The susceptibility does not vary much throughout the whole 
temperature range (Fig. 3 a, b). AC-cleaning was just successful for the 
“soft” components of the magnetization (fig. 3 c, d); the whole collection 
was cleaned at three temperature steps. 
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The NRM-intensities vary from 5 10“^- 10“^ A/m. 

In the case of very low NRM-intensities, a successful isolation of the 
true ChRM was not always possible since their intensity was within the 
noise level of the cryogenic magnetometer. 


Magnetostratigraphic results 

The whole collection of samples were cleaned in three thermal 
demagnetization steps, at 200, 300 and 400 °C. Higher cleaning 
temperatures were chosen, where hematite was clearly identified as main 
carrier of the ChRM. Using all these results a “best fit” column was 
produced (Fig. 4). It clearly shows two polarities; normal at the top and 
the bottom of the section and reversed in the middle part. Since the K/T- 
boundary was defined by geochemical, mineralogical and biostrati- 
graphical methods (Preisinger et al. 1990) the reversed zone is identified 
as G' zone (chrone 29 R), top one as H (chrone 29 N) and the bottom one 
as F 3 '^ zone (chrone 30 N). In some levels a higher scatter can be seen, 
which is most probably due to the turbiditic origin of the sediments. Just 
between 2 and 3.5 m in the Maastrichtian the declination and inclination 
is changing synchronously that a poleexcursion seems to be very 
probable. Any other “finger prints” for the G" zone are missing. 

For the modelling of polarity transition, the interval from 11.8 to 
13.0 m in the Maastrichtian was chosen. As mentioned above, this part of 
the profile consists mainly of limestones and marls. The rocks are very 
suitable for drilling. The cores were drilled with a 30 % overlap of each 
core (Fig. 2). 

Calculating with a sedimentation rate of 1.7 cm/1000 yr. no time 
interval larger then 200-300 yr could have been missed. 


The data underlying the analysis 

Besides the data-set from Gams we used for our modelling two data¬ 
sets from Gubbio (Italy), which were sampled and magnetostrati- 
graphically analysed by Rocchia et al., (1990) Fig. 5, as well as Lowrie and 
Alvarez (1977) Fig. 6. The inclination and declination values of these 
studies are shown in Fig. 5-6. We can see that the F 3 '^-G" polarity reversal 
is normal to reverse transition. 

We worked on both the first and second set of data from Gubbio 
leaving them unchanged, using them without smoothing out anything in 
advance. However, the data by Rocchia et al. show surprisingly small 
scatter, indicating that the data has already been filtered. 

The data of Gams Fig. 7, show a higher scatter obviously due to 
sedimentary processes as mentioned above. 

Non of the possible methods of filtering with respect to the model 
applied leads to useful results. Therefore we left the data of Gams without 
any filtering. 
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Fig. 4: C/T-boundary profile showine the characteristic remanent magnetization with 
declination, inclination ana intensity, susceptibility and polarity. 
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Fig. 6: Inclinadi 
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The model 

Finally to be able to present the data properly, and to develop a 
feeling for the Gaussian coefficients of the transition, we used a flooding 
model, similar to that of Ken Hoffman (1979) and (1981), which is 
discussed briefly. 

Within the core, on a certain point of the lateral surface of a cylinder 
covered with dipoles perturbation is generated. This perturbation spreads 
in both directions at linear velocity, both vertically and azimuthally, 
resulting in a reversal of polarity. Depending on the location of the source 
of the perturbation quadrupolar- and octupolar - dominated transition 
fields can be generated. The variable parameters of the model are the 
dimensions of the cylinder, the location of the source of the perturbation, 
the number of the poles and the vertical step width. The relation between 
the number of poles and the vertical step width can be compared with the 
relation of veolocities (v) which results in a specific value (v-1). 


Results 

Our analysis includes the original, unfiltered data (in general) from 
the range where the inclination finally changes its sign. As a point of 
reference for the model we used the transition of the observed VGP-path 
through the equator. 
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Fig. 8: Model simulation (octupole; perturbation 350° E) for the observed VGP-path from the 
data-set of Lowrie and Alvarez (site: Gubbio; 43.2° E, 12.3° N). 
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It is clearly recognizable that the VGP-path and the site are lying 
within 90° angular distance and can therefore be called nearsided. The fact 
of a nearsided N-R reversal causes a quadrupole transition field, which 
starts in the northern hemisphere or an octupole transition field, which 
starts at high latitudes (Fuller et al. 1979). 

Because of the lack of data from the southern hemisphere, we were 
not able to distinguish between quadrupole and octupole dominated 
fields. Therefore the model VGP-path are nearly the same for octupole 
and quadrupole simulated transitions, but depend in our modelling from 
a different source of perturbation. 

For our model we chose a radius r = 0.13 re (Earth Radii), and the 
length of the cylinder L = 0,8 re. We took the same dimensions of the 
cylinder as K. Hoffman at discribing the Brunhes-Matuyama polarity 
reversal. 

By appropriately choosing the location of the perturbations, we got 
the results shown in Fig. 8-10. - All diagrams show only that part of the 
VGP path in which the transition took place. We achieved satisfactory 
results as far as the octupole-dominated transition is concerned, at a 
perturbation source at 350° east. 

For the quadrupole-dominated transition we chose an perturbation 
source at 310° 



Fig. 9: Model simulation (quadrupole; perturbation 310° E) for the observed VGP-path from 
the data-set of Rocchia et al. (site: Gubbio, 43.2° E, 12.3° N). 
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To achieve a more flat VGP-path we chose the change of the 
dimensions of the cylinder. 

In this case a problem arises, the model VGP-path always starts at 
180° east and ends at 0° so that it is nearly impossible to get a good fit with 
the observed values. In addition a thicker cylinder increases the nonzonal 
coefficients. 

Drawing conclusions as far as the original data is concerned would be 
rather dubious, because the intersection angle of the observed VGP-path 
to the meridians are not supposed to be exclusively the consequence of 
nonzonal non dipole components. 

The Gaussian coefficients of the different transitions are shown in 
Fig. 11-12. For the octupole dominated case, as one can see, the dipole 
(gj°) decreases through zero and the controlling terms are, the octupole 
(g 3 °), the non axisymmetric quadrupole ([g 2 ^]^+[h 2 ^]^)^^^ and the higher 
order term ([g4^]^+[h4^]^)^^^ (see K. Hoffman 1981). 

In case of the quadrupole dominated transition, both zonal 
dipole and octupole terms go to zero while the zonal quadrupole (g 2 °), 
the non zonal quadrupole ([g 2 ^]^+[h 2 ^]^)‘^^ and higher order terms 
are still existent (Fig. 12). 
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Fig. 11: Schematic description of the harmonic content during a octupole-dominated 
transition, generated from the model. 
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Discussion 

Since numerous current investigations on the K/T boundarys will 
lead to large number of reliable data sets of the transition G“ to F % in the 
present paper the first attempt was made based on three available data sets. 
As mentioned above these datasets are coming from Gubbio and Gams 
and have the disadvantage, as far as the modelling is concerned that they 
are close to each other geographically. The currently investigated sections 
will help to overcome this problem since one needs at least 3 data points 
over an angular distance of 180° in longitude for a good estimation. 

Another problem during the modelling is the choice of the cylinder 
dimensions. With regard to the assumption that the cylinder lies within 
the core one gets different possibilities in choosing the dimensions of this 
cylinder. But in this case the goal to determine the distribution of the 
Gauss coefficients with respect to the inclination of the VGP-path to the 
equator was not reached. 

Although all transitions of the observed VGP-path show more or 
less the same sense of inclination to the equator - the three data-sets are 
not enough for such an interpretation. 

The loop in the data-sets of Gams could be due to an excursion before 
the actual reversing process starts, even though there is no evidence in the 
other datasets. The reason for this controversy could be caused by the 
sample distance, in the case of Gams the sampling density was unusually 
high so that much more detail can be expected. Anyway the model 
diagram at the beginning of the reversal shows a close similarity to the data 
previous to the uncontrolled behaviour. This fit is the result of no 
filtering. Any common filtering suppresses the higher inclination values 
and one gets a VGP-path like in Fig. 9. 

Future investigations with the same fidelity as Gams should prove 
whether the loop is a significant fingerprint for this reversal. 

Suggestions that many of the VGP-paths from 100.000 yr to 10 Myr 
are clustered within two meridians have reopened the debate whether the 
transition field is dipolar or non dipolar (Bogue, 1991; Laj et al., 1991). 
However the VGP-path observed from the K/T boundary does not lie 
within these meridional bands. Recently, studies have been made in 
describing the transitional field geometry by using primary (antisymme¬ 
tric about the equator) and secondary (symmetric about the equator) field 
families (McFadden et al., 1988; Hoffman, 1991). In this case there is no 
simple distinction between dipole and non-dipole geometries; each 
family contain both. Maybe this studies will lead to a better description 
of the harmonic content and find out the dominant Gaussian coefficient 
during reversal. The harmonic content of the secondary field family could 
be similar to the distribution of the quadrupole simulated transition. 
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